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Abstract. - We measure stretched exponential behavior, exp (—(t/ro)'^), over many decades 
in a one- dimensional array of coupled chaotic electronic elements just above a crisis-induced 
intermittency transition. There is strong spatial heterogeneity and individual sites display a 
, dynamics ranging from near power law (/3 = 0) to near exponential (/3 = 1) while the global 

' dynamics, given by a spatial average, remains stretched exponential. These results can be 

fSJ ' reproduced quantitatively with a one-dimensional coupled-map lattice and thus appear to be 

. system independent. In this model, local stretched exponential dynamics is achieved without 

frozen disorder and is a fundamental property of the coupled system. The heterogeneity of the 
experimental system can be reproduced by introducing quenched disorder in the model. This 
suggests that the stretched exponential dynamics can arise as a purely chaotic phenomenon. 

^ \ Stretched exponential relaxation is almost ubiquitous in complex systems. It is found 

in glasses Q, spin glasses ||^, polymers, high-dimensional cellular automata random net- 
works m coupled chaotic oscillators [^|^ and others 0. For many of these problems, however, 
' the microscopic origin of this dynamics is difficult to establish both experimentally and theo- 
^ I retically. Understanding the underlying causes for stretched exponential relaxation remains, 
therefore, a central goal in the study of the dynamics of complex systems [^-[ict|, as witnessed 
by the flurry of theoretical, experimental and numerical results that have appeared on the 
subject in the last few years ||l0|"|l3|. 

In the last 20 years, many physical processes have been understood in terms of nonlinear 
dynamics and chaos, these concepts often going beyond the reach of standard equilibrium 
statistical physics. Here, we present experimental results showing a stretched exponential 
dynamics, exp [— (^/''o)'^] i ^ simple one-dimensional chaotic system comprised of a chain 
of chaotic electronic oscillators. In particular, we find that in the vicinity of the transition to 
spatiotemporal chaos, individual site dynamics shows this behavior over six orders of magni- 
tude and that the /3's associated with individual site range from 0.8 (nearly exponential) to 0.2 
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Fig. 1 - Schematic diagram of a single diode-resonator and its coupling to its nearest neighbors in a 
one-dimensional arrangement. 

(nearly power law) depending on the activity of the site. Furthermore, the global dynamics 
of the networks also shows a stretched-exponential dynamics with a f3 having an intermediate 
value. These results can be reproduced using a ID coupled-map model which also displays 
a transition to spatio-temporal chaos. Close to that transition, the model shows uniform 
stretched exponential dynamics. We recover the heterogeneous dynamics of the experiment 
by introducing quenched disorder in the model. These results suggest that the stretched 
exponential dynamics can arise as a purely chaotic phenomenon. 

We study the behavior of a one-dimensional chain of 256 coupled diode-resonators with 
periodic-boundary conditions. Figure |l] shows the schematic diagram of the basic blocks 
forming the one-dimensional array. A diode resonator is comprised of the series combination 
of an inductor (30 mHy) and a rectifier type pn junction diode (1N1004) and is driven by a 
sinusoidal source (100 kHz); it follows the period-doubhng route to chaos as the ac drive is 
increased [|l6|. There are fairly large differences in the individual resonators; for example, the 
drive voltage for a given feature, such as the period-3 fixed point, varies by as much as 20 
%. Diffusive coupling between the elements is provided by resistors as shown and periodic 
boundary conditions are used. This one-dimensional setup has been previously studied with 
unidirectional coupling and under the conditions of stochastic resonance . 

As the drive voltage is increased, the coupled system shows a kink-forming route to spa- 
tiotemporal chaos. At low drive the elements period double nearly synchronously with no 
or very little current in the coupling resistors. Increasing the drive, the system reaches a 
two-band chaotic state in which all elements are in the same band. Following this, coarse 
grained period-2 phase kinks (two-band structures in which neighboring elements are in op- 
posite bands) begin to appear A kink tends to make the system less chaotic since there 
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is now appreciable power dissipated in the coupling resistors and the available phase space 
is reduced. With increasing drive voltage, the number of kinks increases until a regular or a 
near regular kink lattice is formed with a lattice constant function of the coupling resistors. 
At this point each element displays two-band chaos with the phase alternating around the 
chain. Further increasing the drive, the two chaotic bands merge into one, in an example 



of an attractor-merging crisis in a high-dimensional spatiotemporal system |19 . The critical 



drive voltage Vc is where the two attractors just touch. Well above Vc the system displays 
full spatiotemporal chaos. Just above Vc regions of the system still spend long periods of 
time trapped in the old two-band attractor before switching to the new single attractor for 
a time and then back again. In a single chaotic oscillator the intermittent switching between 
two attractors is an example of crisis-induced intermittency and the trap time displays ex- 
ponentially distributed switching times [|9|. As discussed below, crisis-induced intermittency 
also exists in coupled system but the trap time at that point is now dominated by stretched 
exponentials. 

Figure ^ shows a 4096-cycle time series for the 256 elements. The two phases of the 
two-band attractor are represented as black and white, and only every other time step is 
displayed. Note the creation, annihilation and diffusion of the boundaries of the two-band 
attractor. With the coupling resistors used, 150 kil, the coarse-grained spatial periodicity of 
the lattice is about 4 sites. 

In Fig. ^, we plot the distribution of trap times, i.e., how long a site remains in the 
two-band phase. This distribution is measured on a single site over 10^ cycles, or 10 hours 
of experimental measurements. This distribution is very well approximated by a stretched 
exponential with an exponent (3 — 0.40 ± 0.05 and tq = 27 ± 10 drive cycles. In glasses 
and other systems, this behavior is associated with relaxation processes. These processes are 
described by a time autocorrelation function averaged over the entire system. As shown in the 
inset of Fig. H, our system also shows a global stretched-exponential dynamics; the flattening 
of the auto-correlation function at longer time is due to the limited statistics |2^ . 

We can address directly the question of heterogeneities in our setup. Figure | shows the 
activity associated with each site. Activity is defined as the number of times a site changes 
state normalized by the number of time steps. The activity is a measure of the local dynamics 
of the sites: very active sites typically show a dynamics close to that of the chaotic regime 
while the slow sites display a trap-time distribution more in line with the frozen regime. 
Mostly determined by the value of the local chaotic threshold, the variations in activity are 
essentially uncorrelated spatially as can be seen in Fig. IJ. Spatial correlations, measured as 
a two-point function on the lattice, show an exponential decay, exp(— r/rg), with tq — 6 sites, 
in the intermittent phase, close to the basic spatial periodicity just below crisis. 

The level of activity for a given site is also reflected in the overall dynamics of this site. 
Figure ^ shows the trap-time distribution for three sites with different level of activity. All 
three sites are consistent with stretched exponential with exponents /3 = 0.8 for the "hotter" 
site, f3 = 0.4 and P — 0.3, for the "cooler" sites. Faster sites can be equated with an effective 
higher driving voltage and a dynamics moving towards a pure exponential (/? = 1); slower 
sites, on the other hand, display stretched exponential with lower exponent, and can even have 
a power-law correlation ((3 — 0), depending on the external driving. For a big enough system, 
we expect to find the whole range of /3 at any driving in this intermittent phase and see the 
distribution of exponent change with driving, with the value of the global (3 somewhere in 
between the two extremes contrary to what was proposed recently by Richert and Richert jist . 

As is often the case for glasses, where the exponent (3 varies as a function of temperature, 
P here changes with the driving voltage. At low external driving, f3 is zero or nearly zero 
and the system is either totally frozen or showing a power-law distribution of trap times. As 
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Fig. 2 - Time series of the 256 diodes in the intermittent state, just above Vc- Since the system 
displays a fundamental period- two oscillation in time below Vc, we show the state of the diodes at 
even time-steps in a binary representation. 



the driving is increased, P also increases relatively smoothly to f3 — 1.0, for the single site 
chosen, and to about 0.5, on average, before the full system become completely chaotic. The 
relaxation time scale, defined by tq, follows f3 rather closely, roughly doubling in the driving 
range of 3 to 3.2. The overall time scale associated with the traps is therefore of the order of 
100 drive cycles, with trap times reaching many times this value during the measurement. 

Surprising results are obtained in the simulation of a simple one-dimensional coupled-map 
lattice. This model follows Kaneko |M and Johnson et al. U% and is composed of a chaotic 
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Fig. 3 - Distribution of trap times measured over a single site over a 10-hour period by interval of 
1 ms. The raw data are shown to fit to a stretched exponential with f} — 0.4. The apparent tail at 
longer time is due to the discrete nature of the measurement and could be redressed by more statistics. 
Inset: Time auto-correlation function averaged over all 256 sites for a run of of 2048 timesteps of 64 
drive cycles each. The data fits a stretched exponential with /3 — 0.25. 

map, f{x) = 1 — ra;^, coupled to its nearest neighbors in a one-dimensional chain 

Mt + At) = (1 - a),f,{t) + I + ; (1) 

a is the strength of the coupling between the logistic maps. The model, like the experimental 
set-up, shows an attractor-merging crisis with crisis-induced intermittency in the vicinity of 
the transition to spatiotemporal chaos. Below the transition each site is in a coarse-grained 
period-two state. The system also shows spatial periodicity with alternating groups of sites 
in a coarse-grade period-2 as described above for the experimental set-up. The size of these 
groups depends on a. We choose a — 0.25 for a spatial wavelength of 4 sites, similar to 
the experiment. Just above the crisis the system gets trapped in the two-band attractor in 
local regions for a while before going chaotic again. While for an uncoupled system the times 
between switches have a long-time exponential distribution [^9|, Fig. |^ shows that the trap- 
time distribution for this model follows a stretched exponential with (3 — 0.70 ± 0.05 over the 
full distribution, with all sites displaying the same distribution. The stretched exponential is 
therefore not caused by disorder but is a fundamental property of the coupled system. As can 
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Fig. 4 - Activity (fraction of the time steps where a site changes phase) as a function of site number 
for the 256 coupled diode-resonators. 



be seen in Fig. ^ traps show a spatial extension that appears to stabilize them, biasing the 
distribution towards longer time, observed in single oscillators. 

Disorder is necessary to induce the strong site to site variation seen in the experimental 
setup; it can be introduced in the model by selecting a site-varying ri taken from a flat 
distribution of (5, Vi G [r — 5^r + S\. As in the experiment, the global dynamics is also 
well-described by a stretched exponential even though the various sites show a broad range 
in their relaxation dynamics. A fuller discussion of the simulation results will be presented 
elsewhere | |2^ . 

In conclusion, we have found stretched exponential dynamics in a one-dimensional chain 
of chaotic diode resonators near the transition to spatiotemporal chaos. While individual 
sites display a wide range of /3's, the global time autocorrelation function also shows this 
behavior |^^. As can be seen in Fig. || traps can be continuously created and distroyed. They 
are stabilized by their spatial extension, which biases the trap-time distribution towards longer 
times. Based on a simple one-dimensional coupled map model, we show that this behavior is 
not related to disorder. These results indicate that the stretched exponential dynamics can 
be purely chaotic in origin. 

* * * 

We acknowledge partial support from the NSF under grant number DMR-9805848 and the 
Natural and Engineering Research Council of Canada (NM) and the ONR (ERH). We thank 
Don Carter for the device used to measure trap times. NM is a Cottrell Scholar of Research 
Corporation. 



E. R. Hunt, P. M. Gade and Normand Mousseau: Stretched exponential dynamics 7 




50 100 150 200 

t 

Fig. 5 - Trap-time distribution for 3 diode-resonators with different activity, as defined in Fig. ^ 
The main figure shows a log-normal plot of the distribution for all sites. In inset, we show a log-log 
vs. log plot for the three sites; the curves were translated vertically for better display. From left to 
write, in the main panel (and bottom to top in the inset), the corresponding /3 is, respectively, 0.8, 
0.4 and 0.3. 
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